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I. INTRODUCTION 


Although studied intensively for more than 20 years ( e.g ., Giacconi et al. 1962), the 
origin of the diffuse X-ray background radiation (hereafter, XRB) above ~1 keV remains 
controversial. The most popular explanations are thermal bremsstrahlung from a hot in- 
tergalactic medium, with density perhaps of order the closure density, or the superposition 
of a large number of discrete X-ray sources [e.g.. Setti and Woltjer 1979). Strong support 
for the latter explanation has come from observations by the Einstein Observatory. The 
most sensitive Einstein X -ray images (the ‘"Deep Survey” fields) resolve, by direct integra- 
tion of the X-ray logN-logS curve, 26 ± 11% of the XRB into discrete sources (Giacconi et 
al. 19796). Unfortunately, the limiting sensitivity of these deepest of all currently available 
X-ray exposures leaves ~ 75% of the XRB unresolved. 

Einstein observations also strongly confirmed that QSOs were a prime candidate class 
for that discrete source contribution (e.g., Tananbaum et al. 1979). Thus a second approach 
to estimate the discrete source contribution to the XRB has been to use the QSO optical 
logN-logS curve, and a relation (e.g., the mean luminosity ratio < L 7 lL opt >) which 
allows one to predict QSO X-ray emission from that in the optical (e.g., Zamorani et al. 
1981). Although this latter approach is much less direct, potentially it can be used to 
estimate the the contribution of all QSOs, even very faint ones, to the XRB. In order to 
use this optical logN-logS approach, we need to know the X-ray properties of “typical" 
QSOs. By “typical” , we mean QSOs of apparent optical magnitude (and hence, X-ray flux 
and surface density on the sky) potentially capable of contributing substantially to the 
XRB. Unfortunately, the overwhelming majority of previous observations of QSOs with 
Einstein fall into two categories, neither of which are “typical 5 '. 

The first such category is the targeted observations of heterogeneous collections of 
previously known QSOs, many of which are famous for being atypically bright in the 
optical and/or radio (e.g., Zamorani el al. 1981: Ku. Helfand. and Lucy 1980). Although 
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complicating their direct use in the optical log N-log S approach, the heterogeneity (in 2 , B, 
L opt . Lradi selection, etc.) of such collections has led to the realization that < L x / L opt > 
is not constant for all types of QSOs. In addition to a significant dependence on radio 
properties (Tananbaum et al. 1983), < £ 7 /L opt > has been found to depend on L op; 
and/or 2 such that higher redshift and/or higher optical luminosity QSOs have smaller 
< L x / L opt > ratios (Avni and Tananbaum 1982). An issue of much current, controversy, 
and critically important for proper use of the optical log N-log S approach, is which of 
these dependences (L opl or 2 ) is primary (see § 5). Because of such dependences on L opt 
and/or 2 , ensemble distribution parameters such as < L x /L op i > are, in fact, directly 
interpretable quantities only for a sample of QSOs which includes a narrow range of values 
of L opt and/or 2 : such a narrow sample might be used to emprically test various models 
(suggested by the broad samples) for the dependence of < Fj/.Z. op , > on L opt and/or 2 . 

A second category of Einstein observed QSOs are the X-ray selected samples (e.g.. 
Margon, Downes, and Chanan 1985: Gioia el al. 1984, and references therein). However, 
these X-ray selected samples can yield only limited information about typical QSOs; at the 
risk of stating the obvious, X-ray selection tends to preferentially select the X-ray bright 
objects. Further, the dependence of < L T jL opi > on L opt and/or 2 results in a mean 
redshift for X-ray selected samples of only 2 < 0.6. and such low redshift QSOs are too 
rare to contribute substantially to the XRB. 

More recently, however, several X-ray studies of homogeneous optically selected samples 
of QSOs have been discussed by Marshall et al. (1984), Anderson (1985; hereafter ‘“A85"). 
Kriss and Canizares (1985). Tananbaum et al. (1986), Avni and Tananbaum (1986), and 
Schmidt and Green (19S6). These last three papers report on X-ray observations of the 
"XBQS”. an Einstein observed subset of ~ 60 QSOs from the Bright Quasar Sample (BQS) 
of Schmidt and Green (1983). The XBQS is very important because of its unique position 
in the L opt . 2 plane, as well as the high X-ray detection fraction (> 5091). However, the 
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B < 16 QSOs of the sort studied in the XBQS are themselves much too rare to contribute 
substantially to the XRB. Because the BQS is a UV-excess selected sample, it is complete 
only for z < 2.2. Kriss and Canizares (1985) obtained deep Einstein images of two fields 
in which Hoag and Smith (1977) and Sramek and Weedman (1978) had previously grism- 
selected 22 QSOs, most at high redshift. Marshall (1983) obtained deep Einstein images 
of the “Braccesi Faint" (hereafter, BF) sample of ~ 35 QSOs (Formiggini et al. 1980). The 
BF is a UV-excess selected sample complete for B < 19.8 and 2 < 2.2. Thus, although 
the literature contains X-ray information on hundreds of QSOs. X-ray data have been 
reported for only a single complete sample of optically selected QSOs - the 35 objects of 
the BF sample - of the sort that might contribute substantially to the XRB, and there 
have been no such studies of large, complete samples of typical high redshift QSOs. As we 
shall discuss below, X-ray observed complete samples of optically selected QSOs can be 
used with very few uncertain extrapolations to estimate the contribution of such QSOs to 
the XRB. 

In summary, previous data on X-ray emission from typical QSOs lacked the following: 
(l) a narrow sample (in L opt — 2 space) to aid in testing empirically the various previously- 
suggested models for the dependence of < L x /L op( > on L opt and/or z: and, (2) a large, 
complete sample at high redshift to complement the BF sample. In order to obtain such 
samples of typical QSOs. we have taken, in an "after the fact” fashion, grism/grens plates of 
~ 17 deg 2 of the sky previously imaged to very sensitive X-ray fiux levels with the Einstein 
Observatory. The archived X-ray data for our study includes 50 Imaging Proportional 
Counter (IPC) fields and five High Resolution Imager (HRI) fields which have among the 
longest integration times of all Einstein images. The new optical data, obtained at various 
4m class telescopes, include at least one “bhie grism/grens’’ plate centered on each of these 
same deep X-ray fields. In this fashion, we have optically selected more than 400 previously 
uncataloged QSO candidates for which the most sensitive X-ray fiux information is already 
extant. In this report, we will concentrate on a high redshift subsample (~ 90 objects) of 
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this larger data base for which grism/grens selection is highly efficient and reliable. 

In § 2, we present details on the optical selection of our high redshift QSOs. Keeping 
in mind the motivations discussed above, we define two optical samples: a “complete" 
sample and a “narrow" sample. Section 3 presents a discussion of the X-ray data for these 
new samples. Because most of our QSOs are not detected in X-rays, it is necessary to 
incorporate information contained in X-ray non-detections as well as detections; two such 
statistical techniques, “stacking" and “survival analysis’", and some cautions about their 
use are presented in § 4 and the Appendix. Finally, in § 5, we use our data to investigate 
the ensemble X-ray properties ( e.g ., contribution to the XRB, dependence of < L x / L opi > 
on L op t and/or z, etc.) of high redshift QSOs. and the QSO population in general. A 
more detailed account of some of the research described herein may be found in A85, and 
preliminary results of this work have been reported by Anderson and Margon (1983, 1984. 
1986). 


II. THE OPTICAL SAMPLES 

The grisms in use at the 4m telescopes of the Kit.t Peak National Observatory and 
Cerro Tololo Interamerican Observatory, and the “grens’" at the Canada-France-Hawaii 
Telescope, are particularly well suited to discovery surveys for strong emission line QSOs 
(Smith 1981, 1982). The fields of view of the grisms are about 34' x 34'. A CFHT grens 
plate covers about 54' x 54'.; however, due to degradation of image quality near the field 
edges, combined with vignetting, and guide and calibration probe obscuration, the grens 
plates have an effective useful area of order 0.5 deg 2 (Weedman 1985; Crampton, Schade. 
and Cowley 1985). These fields of view are well matched to those of the Einstein X-ray 
images (% 1° x l c ) 

Centered on each of the 55 deep Einstein fields, which are scattered disjointly over the 
celestial sphere, we obtained a IIIaF-j-blue grism/grens plate combination as the primary 
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plate for selection of the QSO sample. The “blue grism" has a blaze wavelength of 3550 A 
and linear reciprocal dispersion of 1500 A mm -1 , while the “blue grens" has 1000 A mm -3 
dispersion with a blaze at 4000 A. Our desire to obtain redshifts directly from the plate 
material was a primary motivation for choosing the IIIaF, with its wider wavelength cov- 
erage of 3200A-6900A. rather than the IllaJ emulsion. Wavelengths of strong emission 
features are estimated from the grism/grens plates to an accuracy of ~ 100 A. To extend 
sensitivity at redder wavelengths, we also obtained “red grism" (blaze ~ 6000 A) plates 
on either a IIIaF emulsion, or on a IV- N emulsion (6100-8900 A with a blue-cut Wratten 
29 filter) whenever time permitted (about 20 fields). In order to maintain as uniform a 
sample as possible for the current study, these red grism plates are used only as additional 
confirmation of the candidate as a QSO. or to confirm the redshift, but not for selection. 
Typically, skv-limited exposures are ~ 45 — 60 min on baked IIIaF plates, and reach a 
limiting continuum magnitude of B % 21.0; exposure times are 60-90 min on the IY-N 
plates hypersensitized in silver nitrate, and reach a similar limiting magnitude in V*. 

The plates were scanned visually using a binocular microscope and an “underlay" grid 
to systematize the search. Both emission line QSO candidates (the main thrust of this 
report) and UY excess objects were searched for. In the ^ 17 deg 2 (effective area) of 
sky surveyed, approximately 400 QSO candidates were selected. Each plate was searched 
several times at separate sessions, and the thoroughness of our search may be addressed 
through the following considerations. A typical limiting continuum magnitude for a blue 
grism-filllaF plate taken at a 4m telescope in better than 2" seeing is B 20.5 (Clowes 
1981). Thirty previously known QSOs (of all redshifts) from the catalog of Veron-Cettv 
and Yeron (1985) lie within the area surveyed by our “high quality’’ plates (all grens plates 
and 38 grism fields, all taken in < 2" seeing), and have B < 20.5; only 8 {27%) were not 
successfully recovered in our “double blind" searches of the plates. More importantly, in 
the redshift regime where grism/grens selection is most successful (r > 1.8). all seven of the 



/ 


previously known high redshift QSOs with B < 20.5 were not only recovered in this fashion, 
but also were assigned the proper redshift. The success with which the previously known 
high redshift QSOs with B < 20.5 were recovered provides considerable confidence in the 
thoroughness of our search technique, to 1 magnitude fainter than the completeness limit 
for which we will argue below. Previously catalogued QSOs which are also the targets of 
the Einstein X-ray images are excluded from the optically selected samples to be defined 
below, even if successfully recovered. This is necessary to avoid potential X-ray biases of 
the sort discussed in the introduction. 

In examining our plates, we have found -- 200 low redshift < 1.8) QSO candidates to 
B < 20 which have strong UV excess (UVX), but which lack any obvious emission features. 
We expect that — 100 of these are indeed low redshift QSOs (see arguments in A85). and 
roughly 40 sucli objects arc the likely optical identifications of Einstein serendipitous X- 
ray sources. However, pending follow up spectroscopy and completeness studies, these UV 
excess objects will not be discussed further. 

a) The. Emission Line Candidates 

There are redshift and magnitude regimes where grism/grens selection is particularly 
efficient and the plates themselves can. with reasonable confidence, be used to classify 
the candidates as QSOs, and estimate redshifts. The efficiency with which emission line 
objects are selected from grism/grens plates depends strongly on line strength and the 
seeing, as well as on the continuum magnitude. The dependence on line strength highly 
favors the selection of “Lyman-a" QSOs at ^ > 1.8. where La is redshifted onto the 
wavelength response of the IIIaF emulsion. With the grism/grens. the spectral resolution 
is determined by the seeing. All CFHT plates were taken in < l" seeing (resolution 
70 A), and we were able to obtain a blue grism plate in < 2" seeing for almost all the 
grism fields (resolution better than ~ 150 A). The frequently excellent seeing at the CFHT. 
coupled with the higher dispersion of the grens. allows somewhat weaker-lined objects to 
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be detected than with the grism (Crampton, Schade, and Cowley 1985). 

The non-uniform response of the IIIaF emulsion as function of wavelength is well 
known. In combination with the blue grism, qualitatively the response is as follows (see 
also Fig. 1 of Carswell and Smith 1978): there is a broad maximum near 3500 A with 
nearly smoothly decreasing sensitivity to 4800 A; redward of 4800 A there are peaks in 
the response at 5300 A, 6100 A, and 6800 A. Because these response peaks can sometimes 
be confused with emission lines, we henceforth will concentrate only on QSOs selected 
from blue grism/grens plates which have an obvious strong feature blueward of 4800 A, 
corresponding corresponding to a < 3.0 for Lo QSOs. The additional spectral information 
redward of 4800 A is used in the current application only to look for confirming spectral 
features. 

b) The Complete .Sample "C 

Clowes (1981) has given a quantitative discussion of selection effects operating in grism 
surveys, concentrating especially on the blue-grism plus IIIaF combination. Clowes' models 
predict, for example, a selection curve for the IIIaF+blue grism plates which provides a 
description of how the plate limiting equivalent width of La+N V varies as a function of 
seeing and continuum magnitude. Clowes' models are quite successful at reproducing the 
appropriate data from the Hoag and Smith grism survey (Hoag and Smith 1977; Osmer 
1980). 

The distribution of rest frame La+N V equivalent widths is available for a large sample 
of radio-selected QSOs (Wilkes 1986). For all blue-grism plates taken in < 2" seeing (38 
fields), we have used Clowes’ models to estimate the limiting equivalent width of La + N V 
detectable to an (observed) magnitude of B < 19.5. Then, by comparing these values 
with Wilkes’ equivalent width distribution, we estimate that our sample of La QSOs from 
these 38 plates should be > 80% complete. A radio selected sample is appropriate for this 
comparison because radio selection is the only efficient line-strength-independent technique 
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of finding high-redshift objects. 

Further confirmation that significant incompleteness sets in at B > 19.5 is provided 
in Figure 1, in which we display the apparent, magnitude distribution for the high redshift 
QSO candidates in our sample. Because the surface density of grism selected QSOs per 
unit magnitude is thought to not decrease until magnitudes of order B ~ 20.5 (Gaston 
1983). the paucity of QSOs to the faint side of the peak of the figure is likely attributable 
to the onset of incompleteness. 

Thus, we have established a subsample (hereafter, the “complete sample", or sample 
“C") of our emission line selected QSOs which should constitute a reasonably complete 
(> 80%) sample. Sample C consists of all emission line QSOs with 1.8 < z < 3.0 and 
B < 19.5 selected from our high quality blue grism plates. When only a single line 
is detected on the plate (true for about half of the objects), we have assumed the line 
identification to be Lo; other workers ( e.g .. Osmer 1980) have found that - 80% of such 
strong-lined grism objects indeed turn out to be Lo QSOs. We argue later that this 
assumption is also a conservative approach when interpreting the contribution of QSOs to 
the X-ray background. For uniformity in optical selection biases, we exclude from sample 
C those QSOs selected on other plate-grism/grens combinations. 

The 38 fields for which high quality blue grism-flllaF plates are available cover an 
area of 12.9 deg 2 ; 39 sample C objects were found in these fields, yielding a surface density 
of 3.02 ± 0.48 QSOs deg -2 to B — 19.5, with redshifts 1.8 < z < 3.0. The 39 QSOs of 
sample C are listed in Table 1. The surface density found here is in good agreement with 
that found by other workers: for example, a value of ^ 2.7 was found in the Hoag and 
Smith survey using the same grism /emulsion combination (see discussion in Schmidt and 
Green 1983). The “complete" sample C will be used in § 5 for a direct estimation of the 
contribution of such high redshift QSOs to the XRB. 


c) The Narrow Sample “A 
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The second subsample of objects we define (hereafter, the “narrow” sample, or sample 
“N”), includes 79 emission line QSOs selected from both high quality grens and grism 
plates, which also fall within the following parameter ranges ranges: 1.8 < z < 2.5 and 
31.0 < log /2500 < 32.0, where /2500 (in units of erg s' 1 Hz -1 ) is the monochromatic 
luminosity at 2500 A in the QSO frame. Sample N does not constitute a complete sample, 
but merely represents those QSOs candidates for which the plate material alone is deemed 
adequate to confidently classify the object as a QSO, and which lie in a narrow redshift 
and magnitude (and hence, optical luminosity) range where grism/grens selection is highly 
efficient. 

In addition, we have also included in sample N another ten QSOs from the the Kriss and 
Canizares (1985) study which also lie within the sample N parameter ranges. These latter 
QSOs were discovered with the identical blue grism used in our study. Sample N, then, 
includes a total of 89 objects which occupy a narrow region of L opt — z space. Although 
incomplete, this sample is nevertheless well suited for empirically testing models for the 
dependence of < L x /L opt > on such other parameters as L op1 and 3, i.e.. for studying the 
X-ray conditional luminosity function (see, e.g., Avni and Tananbaum 1982, 1986). 

The QSOs of sample X have means of optical luminosity and redshift of. respectively. 
< log(/ 2 500 ) > — 31.44, and < 2 >= 2.1, and a listing of the objects is given in Table 2. 
There are 30 objects in common between the samples X and C. 

d) Optical Magnitudes and Monochromatic Luminosities at 2500 A 

In order to determine the ratio L x /L opt . we need to estimate /2500 for these QSOs. 
Optical B magnitudes are estimated using image diameter measurements from the Palo- 
mar Observatory Sky Survey O print, and converted to B magnitudes using the approach 
outlined in Hayman, Hazard, and Sanitt (1979). To avoid systematic biases (personal 
or instrumental), we recalibrated the image-diamcter/magnitude relation for our own 
eye/50x microscope combination (A85). Two “photometric” standard fields were used 
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to test the reliability of our recalibrated image diameter-magnitude relation: the first in- 
cludes 47 photoelectrically calibrated standards from a field in M15 (Sandage and Katem 
1977), while the second includes the 53 grism selected QSOs (Sramek and Weedman 1978) 
with magnitudes determined from spectrophotometry and microdensitometry (Vaucher 
and Weedman 1980). For both these photometric test samples, the rms differences be- 
tween our POSS estimates and the photometric ones are < 0.5 mag for B < 20.7. These 
reliability test results are displayed in Figure 2 and clearly show that any systematic er- 
rors should be quite small in the mean (< 0.1 mag). Our estimates are considerably more 
uncertain (rms ~ 1 mag) at B > 20.7 (A85), and therefore in sample N only those objects 
with B < 20.7 are included. 

The POSS B magnitudes are corrected to continuum magnitudes using the data in 
Veron (1983) to estimate the line contributions, A B (line), to the observed magnitude as 
a function of redshift. The magnitudes are further corrected for extinction within the 
Galaxy, using 

AB(ext) = — 0.24csc \b Il \. (l) 

Finally, the corrected B magnitudes are converted to monochromatic luminosities at 2500 A 
in the QSO rest frame' using relations from Schmidt (1968). Schmidt and Green (1983). 
and Marshall (1983). 

We have assumed a power law spectral form of energy index a c = 1.0 for the optical 
continuum, and adopt <? 0 = 0; Hq = 50 kms -1 Mpc _] unless otherwise stated. These values 
are chosen mainly for ease of comparison wfith previous studies of the X-ray properties of 
QSOs ( e.g ., Zamorani et at. 1981). The optical characteristics of samples C and N are 
showm in Tables 1 and 2. Note that we have revised the / 2 500 values for the ten Kriss and 
Canizares (1985) objects to our assumed a 0 = 1.0. 

e) Optica! Coordinates 

Optical coordinates were measured from the POSS prints using several two-axis mea- 
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suring machines. Although accuracies differ among these machines, all optical coordinates 
should be accurate to at least 5". Such accuracy is more than adequate for the IPC fields, 
but for the HRI with its ~ 3" X-ray error circles, some of these optical coordinate uncer- 
tainties are larger than optimum. Optical coordinates for sample N and C QSOs are also 
listed in Tables 1 and 2. Finding charts for these QSOs are available in A85 or from the 
authors. 


III. THE X-RAY DATA 

The archived Einstein X-ray data for our study are IPC and HRI images which provide 
both spatial resolution and flux information for X-rays of energy a few keY. The IPC has 
a field of view of 60' x 60' (but with ~ 38' x 38' inside the window support “ribs"), and 
spatial resolution of ~ 1.5'. The HRI (field of view of 27' x 27') provides higher spatial 
resolution of 2", but with comparatively low quantum efficiency. 

The particular X-ray images chosen consist of high galactic latitude fields (|5 /y | > 20 : ) 
with very long integration times: 50 IPC images with integration times of 10,000-45.000 s, 
and five HRI images with integration times of 40,000-100,000 s. Thus, the fields we use 
overlap in sensitivity with the so-called “Deep Survey” fields, the handful of the very 
deepest Einstein X-ray images (30,000-60,000 s in the IPC, and 10,000-50,000 s in the 
HRI). and are among the most sensitive of any currently available. 

TIi ese X-ray fields were initially imaged in X-rays for reasons not direct ly related to 
the current study, having instead, a wide range of potential X-ray sources as their primary 
targets, e.g.. bright stars, cataclysmic variables, clusters of galaxies, normal galaxies, and 
AGN/QSOs. However, the X-ray information on the target commonly is contained within 
only the central < 2% of the area imaged. In an “after the fact fashion” (i.e., after the 
X-ray data had already been acquired), we obtained grism/grens plates centered on these 
deep Einstein images. The X-ray properties of the optically selected samples C and X are 
investigated here using the “unused" portions of these X-ray images. 
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The approach used here is to be contrasted with programs aimed at the optical iden- 
tification of serendipitous X-ray sources ( e.g ., Margon, Downes, and Chanan 1985, Gioia 
et al. 1984, Kriss and Canizares 1982, Reichert et al. 1982). The latter approach results in 
an X-ray selected sample of objects. On the other hand, the QSOs investigated here are 
first optically selected from the grism/grens plates, and then the archived Einstein images 
are used to derive a sensitive X-ray flux detection, or more frequently, an X-ray flux upper 
limit. Our optically selected grism/grens samples have no a priori X-ray biases. 

a) Source Detection Algorithms in the 1PC 

The IPC (Giacconi et al. 1979a, 1981) has non-negligible response to, as well as limited 
spectral resolution for. X-rays of energy 0.2-0.28 and 0.4-4. 5 keY. In practice, the objects 
of Samples C and N have such low X-ray fluxes that the IPC must be used in broadband 
fashion. Thus, for the present study, IPC count rates and fluxes will generally be estimated 
for the broadband range 0.5-3.0keV; however, several other broadband energy ranges 
are also used to ascertain source existence. The standard IPC production reprocessing 
(Harnden et al. 1984) has been applied to almost all the relevant fields, and includes 
an automated search technique (in several broad energy bands) to locate X-ray sources. 
For the present study, an optically selected QSO is considered to have been detected by 
the automated source finding algorithms provided that the total number of counts in a 
2.4' x 2.4' region equals or exceeds bo above that predicted for the background, assuming 
Poisson statistics. 

We have also performed a second test for source existence with the IPC data, using a 
larger region of radius 3' centered on the optically selected QSO position: the background 
is estimated from a concentric annular “detect frame" extending radially from 3' to 5'. 
and the 0.5-3.0keV band is used. This circular detect region is sufficiently large that: (l) 
virtually all counts in the 0.5-3.0keY band fall within the 3' circular region for a wide 
range of reasonable possible QSO X-ray spectra: (2) even if the X-ray centroid and the 
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optical position differ by ~ 1', most of the X-ray counts will still fall within the region. This 
latter point is particularly important for “stacking” the images (see below) when the QSO 
is not obviously detected in X-rays. The price paid for using the larger 3' detect region is 
the introduction of extra background noise into the determination of source counts. For 
example, a bo detection of a QSO in the 2.4' x 2.4' region would register only as a % 2.6o 
detection in the 3' radius circle. Thus, in addition to the bo automated detection algorithm 
criterion, we also count as possible detections QSOs which show a > 2.6 o detection in the 
larger 3' circular region. 


b) Fluxes in the IPC 

For the reasons given above, broadband IPC count rates in the 0.5-3.0keY band are 
also determined using the larger 3' circular detect region. Again, this introduces con- 
siderable additional background noise, and a few of the weaker sources (even though bo 
detections in the automated source finding algorithms) have rather large (to ~ 50%) un- 
certainties in their count rates (and hence fluxes). For objects not detected in X-rays, we 
derive a 2.6 o upper limit appropriate to the 3' radius region (corresponding to a bo upper 
limit in the automated detection algorithms). 

The IPC broadband flux F 0 b s (0.5-3.0keV) is derived from the observed source count 
rate as described by Tananbaum et al. (1979). We use a recent recalibration of the IPC 
effective area that causes the X-ray fluxes listed here to be systematically ~ 10% higher 
than the values which would have been found using the earlier calibration. In converting a 
broadband count rate into a broadband flux, it is necessary to assume an X-ray spectrum. 
Following the usual practice, ( e.g Zamorani et al. 1981), we assume that the average 
QSO spectrum (in the 0.5-3.0keY range) is a power law witli spectral (energy) index 
a z = 0.5. This spectrum has been traditionally assumed for it is similar to the value 0.4 
found for the spectrum of the XRB in the 3-20keV band (Schwartz 1979); de Zotti cl al. 
(1982) have shown that if the discrete source contribution to the XRB is substantial, the 
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typical contributors must have X-ray spectra within a few tenths of a x = 0.5. Further, 
such a spectrum is also similar to that of 3C273 (Worrall et al. 1979). However, for lower 
luminosity AGN (particularly Seyfert Is) the mean spectrum between 0.7 and 100 keV has 
spectral index a x ~ 0.7 — 0.8 (Mushotzky et al. 1980, Rothschild et al. 1983, Petre et al. 

1984) . Further, recent work (Elvis, Wilkes, and Tananbaum 1985, Elvis and Lawrence 

1985) indicates that the mean value for a subset of the BQS may be much closer to a x — 
1.0 — 1.2, with a wide dispersion about the mean. The XRB intensity is well measured only 
above ss 2keV, and yet the IPC and HRI have an “effective" energy response of 1-2 keV; 
thus, it is common practice to compare Einstein and XRB data in terms of monochromatic 
fluxes at 2keV. The response of the IPC is such that if the actual X-ray spectral slope 
is of order 1.2 rather than 0.5, the inferred monochromatic flux at 2 keV for a QSO, 
and the contribution of QSOs to the XRB at at 2keV, will be overestimated by ~ 40%. 
This unavoidable uncertainty may dominate our analysis, and we stress that we employ 
a z = 0.5 for conservatism and for consistency with previously published results. The 
combined uncertainties of a,, variations in the IPC gain, and corrections for interstellar 
photoelectric opacity imply a total uncertainty in broadband fluxes of ~ 30% (Tananbaum 
et al. 1979). 


c.) The High Resolution Imager (HRI) 

The HRI (Giacconi et al. 1979o. 19S1) also provides broadband (~ 0.1 — 4.5 keV) 
imaging and flux information. The HRI provides high spatial resolution (~ 2"). but at the 
expense of spectral information, field of view, and sensitivity compared with the IPC. In 
the current application, the primary interest is to push to very faint X-ray flux levels, and 
so only a few very long integration HRI fields have been included in our survey. The lack 
of energy resolution in the HRI introduces added uncertainty into the conversion of count 
rates into monochromatic fluxes. 


Similar analysis methods to those discussed above for the IPC were applied to the 
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HRI images. Sources found by the automated source finding algorithms are checked for 
positional coincidence of < 5" (most of this uncertainty is due to the optical position 
uncertainty) with the sample C and N QSOs, and smoothed X-ray contour maps were 
made. Finally, we followed the approach outlined by Marshall (1983) as a further check. 
At the location of each optically selected QSO, we used a circular region (6.8" radius for 
source < 7’ off axis, and 13.5" for sources > 10' off-axis) to establish source reality; proper 
allowance was made for optical position uncertainties, and the background was estimated 
from a nearby 4' x 4' box which contains no obvious X-ray sources. A detection is again 
defined as bo above background in this circular region. As for the IPC images, HRI count 
rates and fluxes are estimated using a larger region (in this case an 18" radius circle), again 
implying that some genuine sources have rather large uncertainties in their count rates and 
fluxes. 


d) X-Ray Fluxes , Luminosities, ur/d Luminosity Ratios 
for Grism/Grens Selected QSOs 

Broadband X-ray fluxes (corrected for absorption in the Galaxy) are converted into 
monochromatic 2 keV luminosities, / 2 keV assuming a r = 0.5 and using the relations 
given in Zamorani et al. (1981). It is useful to compare the ratio of X-ray and optical 
monochromatic luminosities; this ratio is independent of uncertainties in the cosmological 
parameters Hq and q 0 . and allows for easy comparison of the relative X-ray emissivities of 
QSOs of a variety of optical luminosities. We follow the usual practice ( e.g ., Tananbaum 
et al. 1979), and parameterize this ratio through the quantity a oz . which is the slope of a 
hypothetical power law connecting optical (2500 A) and X-ray (2keY) wavebands in the 
QSO frame, 


l°g (l 2 kev/^250o) 
log^kevM&Oo) 


~ (1°S [h kev//250o)]/2.6O5. 


( 2 ) 


Presented in Table 3 are the X-ray data for sample C and X QSOs. Even at these 



17 


most sensitive of all currently available X-ray sensitivities, only ~ 27% (21/78 for sample 
N, and 10/37 for sample C) of the grism/grens selected objects are detected in X-rays. 
Note that a few sample C and N QSOs from § 2 do not lie within the X-ray image and are 
omitted from the table. The position angles of the X-ray image and the grism/grens plate 
are not identical, so that there are non-overlap regions in the corners, for example. 

At the very sensitive X-ray and optical flux levels of our data, the chance of a positional 
coincidence between a grism/grens selected QSO and a physically unrelated X-ray source is 
non-negligible. Using the X-ray logN-logS curve from the Deep Survey study of Giacconi 
et al. (19796), we predict that to a typical limiting sensitivity of 5 x 10~ ] ^ erg cm -2 s _1 
(0. 5-3.0 keV), there is of order one such chance coincidence within samples C and N. 

Recalling that 43% of the comparably optically faint BF QSOs were detected in X- 
ravs (Marshall 1983), the low rate of positive flux detections reported here gives the first 
qualitative indications that our high redshift objects might have lower ratios of X-ray to 
optical emissivities than their lower redshift/lower optical luminosity counterparts. One 
way to parameterize the mean X-ray to optical luminosity ratio of a sample of QSOs is 
the parameter aS? (Tananbaum el al. 1979). By analogy with the definition of a oz above, 
this parameter is defined as 

Qoi = -[log < 1 2 keV'/^2500 >]/2.605. (3) 

where < hkcx/hboo > is the mean of the ^Vex/hboo ratio for the sample. Note that aS? 
is not the same as the mean a ox , < a ox >, which we also employ elsewhere. In order 
to determine these mean parameters for samples C and N, and quantitatively compare 
them with other samples, it is essential to use a statistical approach which incorporates 
information contained in the X-ray non-detections. 

IV. TWO STATISTICAL APPROACHES FOR INCORPORATING 


INFORMATION IN X-RAY NON-DETECTIONS 



IS 


In this section, we briefly discuss two techniques to deal with “censored" data, i.t., data 
with some upper or lower limits. The first approach, “survival analysis", has been widely 
used in X-ray studies of QSOs. Survival analysis is powerful in the sense that it allows 
an estimation not only of single ensemble parameters like aS?, but also of, for example, 
the distribution function of a oz values. However, we show in the Appendix that in certain 
circumstances and in the presence of a substantial non-random censoring, an assumption 
commonly made in survival analysis can lead to misleading results. The second technique, 
“stacking”, has not had any widespread usage in such studies of QSOs. This alternate 
approach serves as a useful check on the results of survival analysis for the estimation of 
some ensemble X-ray parameters; however, it does not serve as an estimator of distribution 
functions. 


a) Survival Analysis 

Survival analysis was introduced into widespread astronomical use for precisely the 
current sort of application: Avni et al. (1980) developed a technique (“detections and 
bounds”) for incorporating information contained in lower limits, primarily with the in- 
tention of investigating the distribution function of a oz for the famous QSOs which were 
Einstein targets. Pfleiderer and Krommidas (1982) reformulated this technique for studies 
of extragalactic radio sources, using an information theory approach. More recently, it has 
been realized that the “detections and bounds" method is an independent redevelopment 
of a statistical approach known as survival analysis which was already in common use in 
the fields of epidemiology, actuarial science, and industrial reliability (Feigelson and Nel- 
son 1985, Schmitt 1985, and references therein). These latter references also provide an 
excellent introduction to the field with an emphasis on astronomical applications. 

In the current application, there are objects detected in X-rays for which a oz is known, 
and lower limits on a OI for other objects. A convenient, non-parametric way to estimate 
the distribution of a random variable ( e.g ., q oz ) in the presence of “right" (lower limit) 



19 


censoring has been given by Kaplan and Meier (1958), who provide a maximum likelihood 
derivation. We desire to estimate the distribution function of a ox values for a sample of n 
QSOs. Let S(a 0I ) be the so-called “survivor function”, the probability of finding a QSO 
with a > a ox . Let a, be the the true value of a ox for the ith QSO, and a| im be the 
corresponding detection limit (this will vary from QSO t.o QSO). If a, < the object is 
detected in X-rays; otherwise it is only known that a, > aj' m . Then let a t = min (a,-, aj ,m ) 
be the measurement of (or limit on) a OI for the fth QSO; for simplicity, assume the 
measurements are ordered so that a, < a 2 < ... < a,. ..a„. Then the Kaplan-Meier 
product limit estimator of the survivor function can be shown to be 


s (“~> = n (> ~ f ) • 


( 4 ) 


where n t (a,) is the number of objects in the sample with measurements (detections or lower 
limits) a > a.,. dj(a,) is the number of detections at a r , and the product is over values of t 
such that a, < a oz . If the largest observation is censored, then Eq. 4 shows that S(a 03 .) is 
finite as a 07 approaches infinity; to avoid such an unphysical result, the largest observation 
is commonly treated as uncensored, even when it is censored (e.g., Feigelson and Xelson 
1985). 

The mean value of a OJ is estimated from 


/•oc roc 

< a ox >= / af(a)da = / S(cx)da. 
■I o Jo 


( 5 ) 


n-r 1 

= V s ( a >) K - 

7=1 

(Feigelson and Nelson 19S5). If a „ is censored, 5(oc) > 0. and < a OI > is formally 
infinite; hence, the usual practice is again to treat a„ as a detection even if censored. In 
the Appendix, we show- that, for certain non-random censoring patterns (in particular, in 
the presence of a sensitivity threshold), this common assumption can lead to misleading 
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results. However, this potential pitfall has not significantly compromised published results 
on X-ray emission from QSOs. Analytical estimates of the uncertainties in the mean of 
the sample and the survival function are also available (Feigelson and Nelson 1985 and 
references therein). 


b) Stacking 

Given the sorts of uncertainties noted in the Appendix for survival analysis, it is 
desirable to have a second technique to check the former results. Caillault and Helfand 
(1985) have innovatively applied stacking techniques to X-ray images of the Pleiades in 
order to estimate the mean X-ray luminosity of Pleiades stars. We have applied this 
technique to our samples C and N. 


i) Mean Flux for Complete Samples 

In the X-ray image at the location of each optically selected QSO, the total counts 
(0. 5-3.0 keV) in the 3' radius region, and the background counts in the 3 — 5' annulus 
are determined. These counts are summed [i.e., stacked ) for all QSOs in the sample, and 
the mean flux of the sample is estimated. For a sample of n QSOs with X-ray images of 
comparable flux sensitivity, this gives the usual 77 3 ' 2 increase in sensitivity over that in a 
single X-ray image. 

Let 5, be the X-ray source counts (0. 5-3.0 keV) for the f-th QSO. so that the X-ray 
flux (0.5-3.0keV) for the f-th QSO is F, = C, S, / ■” 1 . where t, is the integration time of 
the appropriate X-ray image, and C, is the conversion factor for transforming broadband 
count rates into broadband fluxes. Let Tj and B , be. respectively, the total (background 
plus source.) counts on the source and background counts for the i-th QSO. Then, it is 
readily shown (see A85 for details) that the mean X-ray flux for the sample of n QSOs is 
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just 


< F >~< C > 


E«.- 


=< c > 


( T-B ) 

t 


(6) 


where T. B . and t are, respectively, the total counts, background counts, and integration 
time in the stacked image. Further, the uncertainty in this mean flux is (again see A85) 
just 


0<F> — < C > 


{T + B ) 1/2 

_____ 


( 7 ) 


In deriving Eq. 7 we have assumed the following: (l) Poisson statistics; (2) weighting by 
the inverse of the variance; (3) the background count rate from image to image to approx- 
imately constant; (4) either that the individual object count rates are small compared to 
the background count rates, or that the count rates for all QSOs are approximately equal; 
and (5) that C t is approximately constant for all QSOs in the sample. 

Thus, although the individual 5,’s cannot be accurately determined (for low X-ray 
fluxes) from the single images, the ensemble average flux can be sensitively measured from 
the stacked image. The ensemble mean X-ray flux is the quantity required for estimating 
the contribution of a complete sample of QSOs to the XRB; for this application, it is largely 
irrelevant how the X-ray flux is distributed among the individual objects. We apply the 
stacking technique to sample C (§ 5) to obtain such an estimate for QSOs to B < 19.5 and 
1.8 < z < 3.0. The possible consequences of flux thresholds on the stacking technique are 
considered in the Appendix. 


ii) Stacking and a|5 for Sample N 

Stacking can also be used for sample N to estimate a£?; and this will again serve as a 
useful check on the results of survival analysis. Let r = / 2 ke \ 7 /^2500 be the monochromatic 
luminosity ratio. Then 

r = (1 + 

72500 


( 8 ) 



where /2keV and /2500 are the monochromatic fluxes in erg cm -2 s _1 Hz -1 in the observer’s 
frame at the respective frequencies. Then for a sample of objects like sample N which 
occupies a narrow region of L opt . z space and which is imaged in X-rays within images of 
comparable sensitivity, it can be shown that (see A85) the ensemble mean ratio is 


<0(1 + <2 >) Q *- a ° £ Si _ < C > (1+ < 2 >) Q *- Q ° [T -B) 

r V f* 1 ~ V' p 1 ! ' ' 

2500 ' » ^2 2500 + 

where T and B have the same meaning as above for the stacked image, but where < C > 
is now a conversion between broadband counts and monochromatic flux. An estimate of 
the uncertainty in < r > is 


°<r> 


< C > (1+ < 2 (T + B ) 1/2 

^ /2500 


( 10 ) 


neglecting dispersion introduced by optical flux uncertainties, and by using mean values 
< z > and < C >. Then, an estimate of aS? is given by = — log(< r >)/2.605. 
Note that the numerator in Eq. 9 is proportional to the total number of X-ray counts 
emitted jointly from all the QSOs in time t. while the denominator is proportional to the 
total number of optical counts received in the same time. In § 5 we will use this stacking 
estimate on sample N to check the results of the survival analysis estimate. 


V. X-RAY EMISSION FROM HIGH REDSHIFT QSOs, AND THE 
CONTRIBUTION OF QSOs TO THE DIFFUSE X-RAY BACKGROUND 

a ) The Contribution of B < 19.5 QSOs to the XHB 
For an optically complete sample with X-ray information available, one can directly 
estimate the contribution of such QSOs to the XRB merely by estimating the mean X-ray 
flux of the sample, and in turn multiplying by the surface density of objects in the complete 
sample (Marshall et at. 1984). The details about how that X-ray flux is distributed among 
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the individual objects and how a ox depends on other parameters is largely irrelevant for 
such an application. The 37 sample C QSOs in Table 3 are used in this fashion. The current 
approach, using a complete sample, is to be strongly contrasted with most optical logN- 
log S estimates of the contribution of QSOs to the XRB. In their Table 2A, Maccacaro, 
Gioia, and Stocke (1984) list ~ 21 parameters or uncertain assumptions required in using 
the usual optical logN-logS approach. 

The approach used here with an optically complete sample requires only about six 
such parameters or assumptions: a complete sample is characterized by the redshift range 
where it is thought to be complete, the completeness magnitude limit, and a surface 
density to that limit, and of course one must assume a value for a 2 (this latter assumption 
is the primary uncertainty). With only a half-dozen or so parameters or assumptions 
involved, this approach is at least competitive with the X-ray luminosity function approach 
emphasized by Maccacaro, Gioia, and Stocke (1984). Faint, optically complete samples 
imaged to faint X-ray flux levels are potentially powerful estimators of the contribution of 
QSOs to the XRB. 


i) Survival Analysis of Sample “C" 

Application of survival analysis to the 37 sample C objects with X-ray information 
available yields the following estimate for the mean monochromatic flux at 
2 keV: < 5 2 j te v >= (1-24 ± 0.10) x 10“ 3 keV cm -2 s -1 keV “ 1 (corresponding to 
< jF( 0.5 — 3.0 keV) >= 3.97 x 10~ 14 erg cm -2 s _1 in the broadband). The fractional 
contribution of such QSOs to the XRB is then, 

rQSO . o 

f(B < 19.5, 1.8 < a < 3.0) = = 0.021 ± 0.004, (11) 

■ / 2keV "*2 keY 


where n is the surface density of the high redshift objects derived for the sample C (3.02 ± 
0.48 deg~ i) 2 ), and = 5. 84 keV cm -2 s _1 ster is the XRB intensity at 2keY (Schwartz 
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1979). The error estimate of Eq. 11 includes both the uncertainty in the mean flux and 
that due to Poisson fluctuations in the surface density of QSOs. 

ii) Contribution to the XRB of Sample C QSOs Using the Stacking Approach 

We have also used stacking to estimate the contribution of sample C QSOs to the XRB. 
In the stacking technique, it is necessary to exclude QSOs for which the detection region 
is strongly contaminated by a nearby but not physically associated X-ray source (objects 
17272 -1- 502. 17274 -I- 503. and 21556 — 302). We also must exclude objects imaged in the 
HRI, since, for instrumental reasons, they have very different (compared with the IPC) 
values for Cj and (7) 4- B,)/tj (§ 4: A85); thus, Obj. 11, Obj. 15, Obj. 16, and 14149 + 251 
are also excluded in the stacking approach. 

Application of the stacking method to the remaining 30 class C QSOs yields for the 
stacked image 7=4348 counts. f?=3821 counts, /=5261 37 s, and thus, S = (T - B) = 527 
counts, 6S = (7 + B) ]l 2 — 90 counts. < 5 >= S/t = 1.00 x 10" 3 counts s'" 1 . The 
mean conversion between broadband count rate and fluxes for these objects is < C > = 
4.34 x 10" 11 erg cm -2 s -1 per count s _1 . which yields 

< 7(0.5 — 3.0 keY) >= 4.34 x 10“ H erg cm" 2 s" 1 (12) 

or, 

< S 2 keV >= (9.4 ± 1.9) x 10" 6 keV cm" 2 s" 1 . (13) 

Included in Eq. 13 are the expected uncertainties due to the assumption of the constancy 
of the the terms C, and (7, + B^jt,. Effectively, this stacked image has a limiting flux of 
2.6 < C > B ] i 2 /t = 1.35 x 10" 14 erg cm" 2 s" 1 (0.5-3.0keY ), comparable to the values 
reached in the most sensitive “Deep Survey' 1 fields of Giacconi el al. (19795). 

Thus, for the fractional contribution of sample C QSOs to the XRB, the stacking 
technique yields: 


f(B < 19.5, 1.8 < z < 3.0) = 0.016 ± 0.004. 


M 
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This value is within la of that found here using survival analysis. Given the very different 
approaches, it is comforting that these two estimates are in such good agreement. The 
fractional contribution found here for sample C QSOs is quite small, and results in part 
from the low surface density of such high redshift objects, but also from the comparatively 
low value of X-ray emissivities from these objects (see discussion below). 

Hi) Contribution of all QSOs with B < 19.5 to the XRB 

Sample C is an excellent complementary sample for combination with the BF sample. 
We have applied survival analysis to the 14 BF QSOs which have B < 19.5 and z < 1.8 
(and with Mb < —23.0). We have not used the stacking technique with the BF objects 
because the actual counts are not published for the non-detections. For this BF subset, 
we find 

< i> 2 keV >= (3.4 ± 1.4) x 10~ 3 keV cm -2 s" 1 keV -1 . (15) 

Then to B < 19.5, the contribution to the XRB from QSOs in the range 0 < z < l.S 
is 

f(B < 19.5,0 < s < 1.8) = 0.16 ± 0.08. (16) 

Thus, even using the higher survival analysis value (Eq. 11) for the high redshift QSOs 
and allowing for ~ 20% incompleteness in sample C. and 11% in the BF sample (Marshall 
1983), we estimate the XRB contribution from QSOs to B < 19.5 and 0 < z < 3.0 to be 

f(B < 19.5.0 < z < 3.0) = 0.20 ± 0.08. (17) 

The surface density of very high redshift QSOs (s > 3.0) is the sxibject of a number of 
recent and ongoing investigations (Koo and Kron 1980, Osmer 1982, Hazard and McMahon 
1985, Schmidt, Schneider, and Gunn 1986). Although the situation is far from settled, it 
seems likely that unless such QSOs have efficiently escaped selection by current methods, 
such objects probably have surface densities to B < 19.5 which (conservatively) are no 
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more than ~ 20% of the surface density of QSOs in the redshift range 1.8 < z < 3.0 for 
B < 19.5. It therefore seems quite unlikely that, a > 3.0. and B < 19.5 QSOs contribute 
more than ~ 0.4% of the XRB. A reasonable estimate for the contribution of all QSOs to 
B < 19.5 to the XRB is thus 


/(£< 19.5) = 0.21 ±0.08. (18) 

As noted earlier, there are a number of reasons to believe that the values estimated 
above for the contribution of (B < 19.5) QSOs to the XRB might be more properly re- 
garded as upper limits : (l) possible contamination of Sample C by lower redshift (hence 
lower L opt , lower a ol ) QSOs; (2) a z may be considerably steeper than 0.5; (3) threshold ef- 
fects in the survival analysis estimates; (4) probable chance positional coincidence between 
the order of one sample C QSO and an unrelated X-ray source: (5) possibility that the BF 
surface density for B < 19.8, 2 < 2.2 QSOs is abnormally high (see Marano, Zamorani, 
and Zitelli 1984 for a discussion of this controversial issue): and, (6) ~ 40% of the total 
X-ray flux from the B < 19.5. z < 1.8. subsample arises in a single object, BF 222. 

Contamination of the sample C by non-QSOs of course would tend to dilute the en- 
semble average X-ray flux derived for the sample. However, this diluted X-ray flux is then 
multiplied by an inflated surface density and these two effects should roughly cancel. 

b) The Distribution of a ox Values for Sample N QSOs 

Avni and Tananbaum (1982, 1986) have empirically found a dependence of < a 0I > 
on L op t and/or 2 ; which dependence is primary is of current controversy and vet is crucial 
to proper use of the optical logN-logS approach. Using a collection of optically selected 
QSOs and Sevferts spanning a large range of L op i and 2 , they find a primary dependence 
on L opf using linear regression techniques appropriate for censored data: 


< a OI >— — 0.0l[r(~) — 0.5] -f- O.llSjlog / 2500 


30.5] + 1.48. 


(19) 
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where r(z) = z/( 1 + z) is the fractional lookback time (for go = 0): and where we have 
corrected this expression for a recent IPC effective area recalibration. Such a dependence 
on L opi can imply substantial (> 50%) contributions to the XRB ( e.g ., A85, Marshall et al. 
1984). and the relationship of Eq. 19 has been given various physical interpretations (Tucker 
1983, Netzer 1985). However, Schmidt and Green (1986), in trying to reconcile data from 
the XBQS with data from the Einstein Medium Sensitivity Survey, favor evolution with 
redshift, such that L z fL opt ~ (1 4- z)~ i) 2 . They find values of order only 10% for the 
contribution of QSOs to the XRB. Selection biases that correlate L opt with 2 have made 
it difficult to confidently disentangle the effects of / 2500 and 2 on < a OJ >. 

Because of such dependences on L op t and/or 2 , ‘"observable" quantities such as < a oz >, 
aS?, and the survivor function for a oz values are immediately interpretable only for a 
sample of QSOs which occupies a narrow region of L opl . 2 space. Sample X occupies 
a sufficiently narrow region of L opt , 2 space that it can be treated as an empirical data 
point for QSOs with optical luminosity and redshift equal to the means of the sample: 
(< log (^ 2500 ) >= 31.44) and (< 2 >= 2.1). 

i) The Survival Analysis of Sample A 

For sample N, the survival function, not merely a single parameter such as < a oz > 
or ctS?, is of interest, and so it will be generally more productive to use survival analysis 
rather than the stacking technique. Using the Kaplan-Meier product limit method of 
survival analysis, we have estimated the distribution of a OJ values for the sample X QSOs 
shown in Figure 3. For sample N the survival analysis estimates yield 

< a OT >= 1.58 ± 0.03 and af z = 1.50 i 0.03. (20) 

These high redshift, high L op t QSOs are indeed comparatively weak in their relative 
X-ray emissivities. For example, the X-ray selected QSOs (e.g.. Margon, Downes, and 
Chanan 1985; Gioia et al. 1984) as well as radio bright QSO (e.g.. Zamorani et al. 1981. 
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Tananbaum et al. 1983) both have a% x ss 1.3; the sample N QSOs have a mean < L x j L opt > 
ratio that is ~ 3 — 4 times smaller than that for the X-ray selected and radio bright QSOs. 

For the mean parameters, < log (/ 2500 ) > = 31.44 and < 2 >= 2.1, appropriate for the 
sample N QSOs, the model of Avni and Tananbaum (Eq. 19) predicts values of < a 0I > = 
1.59, and aS? = 1.47. Both values in Eq. 20 are within lcr of these model predictions. 
Thus, the measurements of < a 07 > and o.§? for sample N strongly confirm the viability 
of the relationship found by Avni and Tananbaum for the dependence of < a ox > being 
(primarily) on L op t . Note that the measurements here have made no assumptions about 
the particular parametric form (or even the existence!) of the dependence of a ox on L opt 

and/or 2 . The measured values of a ox and a„ x are also independent of uncertainties in JT 0 
and q 0 . 

In order to investigate the possibility of redshift evolution of a oz , we have compared 
our a ox data for sample N with that for various subsets of the XBQS. These subsets of 
the XBQS are chosen to span the same range in L opt as that spanned by sample N, and 
the different subsets are appropriate for different assumed values for qo • These XBQS 
subsets also span narrow redshift ranges, but at a very much lower mean redshift (typ- 
ically < 2 0.3 — 0.4 rather than < 2 >~ 2 for sample N). We use the Mantel-Cox 

and Breslow-Wilcoxian tests ( e.g ., Feigelson and Nelson 1985) to test the equality of the 
distribution for sample N with each of the XBQS subsets. The results are shown in Figure 
4. For ease of comparison with the XBQS data, here (and only here), we have assumed 
q 0 = 0.5. The formal statistical tests confirm what is apparent to the eye: there is lit- 
tle evidence for redshift evolution. Both tests support the hypothesis that, sample N and 
the XBQS subset arise from different parent distributions only at < 35% confidence level 
for the cases of <? 0 = 0 and 0.5. The tests performed here are non-parametric, model- 
independent empirical evidence against redshift evolution of a 0T . Thus, our data support 
models in which the primary dependence of a OJ is on L opt rather than redshift. The pri- 
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mary uncertainty remaining in this conclusion is the previously mentioned possibility of 
modest contamination of sample N by low 2 QSOs. 

Since the values for < a ox > and found here for sample N are in good agreement 
with the Avni and Tananbaum model of Eq. 19, we have further used their model to correct 
the observed sample N distribution for the small (but non-zero) spread in L op t- Following 
Avni and Tananbaum, we determine, for each sample N QSO, the value (or limit) on the 
quantity (see Eq. 19) 

ft^ rr = G 0 i - 0.118 (log l 2h 00 - 31.44). (21) 

The survivor function for these “corrected" values of a OT is shown in Figure 5 and, the 
corrected values of < a ox > and aS5 are 1.56 ± 0.03 and 1.50 ± 0.03, respectively. As 
expected, these values are quite similar to those obtained directly from sample X without 
any corrections, as can be seen by comparison of Figures 3 and 5. Again, we use the 
Mantel-Cox and Breslow-Wilcoxian tests to compare the corrected and uncorrected a 0I 
distributions for sample N. The hypothesis that these two distributions actually arise from 
distinct parent distributions is supported only at the < 45% confidence level. 

ii) The Width of the a 0I Distribution for Sample N 

Franceschini, Gioia, and Maccacaro (19S6). Zamorani (1985). and Avni and Tanan- 
baum (1986) have recently emphasized that the intrinsic width of the a OI distribution. 
o(a ox ), about the mean relation of Eq. 19 is an important parameter for understanding 
the X-ray properties of the QSO population. In particular, the large value (~ 0.2) found 
by Avni and Tananbaum for the standard deviation of the a ox residuals about the mean 
dependence on L opi . together with an assumed optical luminosity function, leads to some 
inconsistencies between X-ray and optically selected samples of QSOs. For example, using 
the optical luminosity function and the a OI distribution for optically selected QSOs, one 
can predict the X-ray luminosity function. However, this process with a standard devi- 
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ation of o(a 0z ) « 0.2 found by Avni and Tananbaum (1982) overpredicts by a factor of 
~ 2 — 3 (Franceschini, Gioia, and Maccacaro 1986; Avni and Tananbaum 1986) the number 
of X-ray sources that, should have been found in the X-ray selected Medium Sensitivity 
Survey (Gioia et al. 1984) sample. Indeed, the desire to resolve such a discrepancy was a 
primary motivation for the Schmidt and Green (1986) model of u 0I redshift, evolution. 

Franceschini. Gioia and Maccacaro (1986) suggest that, an intrinsic value of o(a 0Z ) ~ 
0.14 could resolve the apparent discrepancy, and suggest a number of possible external 
effects that might account for the broadening of the intrinsic distribution, e.g., optical and 
X-ray variability. In fact, for sample N, we find an observed value of a(a ox ) ~ 0.14. For 
the corrected distribution, the standard deviation about the mean dependence of Eq. 19 
is o(olo° tt ) = 0.12. 

Superficially, these small standard deviations found for sample X appear quite inter- 
esting. However, in the Appendix we show that survival analysis threshold effects can 
produce an artificial truncation of the distribution which spuriously narrows and skews 
the survival analysis estimated distribution. It is really the width of the low side of the 
a oz distribution which is most relevant; the contribution to the standard deviation for a oz 
less than < a oz > is 0.13 and 0.11, for, respectively, sample N and the corrected sample 
N. Thus, if the true a oz distribution is approximately Gaussian, the standard deviation in 
the actual sample N distribution is again likely to be ~ 0.2, in agreement with that found 
by Avni and Tananbaum. Otherwise, the actual distribution must be quite skewed toward 
low a oz values. 


iii) Stacking Estimate of aSf for Sample N 

We also have employed the stacking technique to estimate for sample N. As with 
sample C, there are a number of objects not suitable for inclusion in the stacking process: 
the appropriate data for the 10 QSOs from Kriss and Canizares are not published: six 
objects are in HRI fields (00151 + 160. 12275+024. 14149+251. 15090-092, Obj. 11. and 
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Obj. 15); five are severely contaminated by strong, nearby, but unrelated X-ray sources 
(03061 + 169, 11143 + 184, 17274+503, 17272+502, 21556-302); and two are excluded be- 
cause they are not imaged to comparable limiting X-ray flux sensitivities (15100-089 and 
17252+499). Application of stacking to the remaining 55 suitable sample N QSOs yields 

aS? = 1.53 ± 0.02. (22) 

This value is within ~ 1 a of the survival analysis estimate for sample N QSOs. 

c) The Contribution of QSOs to XRB 

Using an optical logN-logS curve, and the a 0J dependence on L opt confirmed by our 
data, one can, in principle, now estimate the contribution of all QSOs to the XRB. However, 
this involves the introduction of the lengthy list of uncertain assumptions and parameters 
noted above: note that none of the proceeding results rely on this uncertain approach. 
Nevertheless, we can at least ask what sorts of contributions are consistent with these 
data. 

For this exercise, we use the optical logX-logS curve of Setti (1984), which he derived 
using the following QSO samples: the BQS {B < 16.0, and z < 2.2): the li BFG” sample 
of Brac.cesi, Formiggini, and Gandolfi (1970) (a color selected sample to B < 18.3. and 
2 < 2.2); the BF (B < 19.8, and a < 2.2); the Kron and Chiu (1981) proper motion selected 
sample ( B % 21.1). and the two-color selected sample of Koo and Kron (1982) ( B k 22.6). 
In deriving this curve. Setti applied incompleteness corrections by, for example, also using 
data on high redshift QSOs from slitless spectroscopy surveys, to obtain the following 
curve for QSOs of all absolute magnitudes: 

log N{< B) % 1.58 + 0.91(1? — 20) for B < 20.1 (23) 


log A r (< B) % 1.67 + 0.21(1? - 20) for B > 20.1 
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where N{< B ) is QSO surface density (counts deg -2 ). We have made the (uncertain) 
assumption of a sharp break at B = 20.1; a smoother function is a less conservative 
assumption in this context. Although derived using only counts to B ~ 22.6, this model 
also agrees well with a variety of more recent work (Hamilton 1984, Boyle et al. 1985, 
Borra and Lepage 1986), some of which extends to even fainter magnitudes. 

The combined X-ray intensity from all QSOs with magnitudes in the range B^^ to 
B max . i n units of erg cm -2 s -1 Hz -1 ster -1 , can be expressed (Zamorani et al. 1981) as: 


r QSO _ 
■*2 keV ~ 


2500 

4410 




(1 



10 -0.4S 10 - 19.380 10 — 2.6D50S? 


dN 

JB 


dB. 


(24) 


where z e g is the redshift of a typical QSO contributing to the XRB. Here, we define 
a QSO to have Mb < -23.0. Then B max is approximately set by the constraint that 
M/j(H max , 2 e ff) < -23.0. This latter constraint is equivalent to assuming that all QSOs 
are at the same redshift 2 = .s eff . We take B ^ n = 16.0, since brighter QSOs are so rare 
that their contribution to the XRB will be entirely negligible. 

By analogy with the Avni and Tananbaum model of Eq. 19, we presume that a£j 
depends on L opi such that, 


aSz = aS? (31.44) + 0.118(log / 25 oo — 31.44). (25) 

where aS? (31.44) = 1.50 is the value found here for sample N QSOs with < log(/ 2 5 oo) > = 
31.44 and < z >= 2.1. Then combining Eqs. 23 (extrapolated to as faint as B ~ 26), 24, 
and 25 yields the results in Table 4. 

As noted above, the entries in the Table are calculated assuming all QSOs to be at z e ff. 
Then, in principle, the actual QSO contribution to the XRB is found by multiplying the 
Table entries by the fraction all QSOs, as a function of B, that are actually at £ = z e {\, and 
summing over the rows and columns. Unfortunately, especially at very faint magnitudes, 
these fractions are not well known: approximately, however, the BF sample, sample C, 
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and the very faint sample of Koo (1986), suggest that z e fj < 2.5 to B < 22.6. Thus. 
QSO contributions as large as ~ 70% are consistent with these data. The additional 
contributions of lower luminosity active galaxies and clusters of galaxies might readily 
account for the remaining ~ 30% of the XRB (Fabian 1981, Schmidt and Green 1986). 

Also note from Table 4 that, whatever the absolute fraction of the XRB due to QSOs, 
roughly half the contribution from QSOs is likely to arise from apparent magnitudes brighter 
than B ~ 20. Such QSOs are the sort directly studied in X-rays in the BF sample (Marshall 
et al. 1984) f or z < 2.2, and in this report for z > 1.8. Qualitatively, we also infer the 
following from Table 4. If QSOs are the dominant contributors to the XRB, the typical 
contributor is probably a low to moderate optical luminosity QSO (apparent magnitude 
near the turnover in surface density counts at £ ~ 20) and at a redshift of z ~ 1 — 2. High 
redshift QSOs apparently are too scarce, and have too low L z /L opl ratios, to contribute 
substantially, while QSOs at z < 0.5 are simply too rare to make a substantial contribution. 

VI. SUMMARY 


We have presented sensitive X-ray information for ~ 90 previously uncataloged QSOs 
in the redshift range 1.8 < z < 3. Even with the longest existing Einstein Observatory 
X-ray exposures, only 25% of these objects are positively detected in X-rays. Among our 
primary conclusions are the following. 

(a) For B < 19.5 and 1.8 < z < 3.0, our sample C data imply that the fractional contri- 
bution of such QSOs to the XRB is less than 3%. 

(b) Combining the results for sample C with X-ray data for the BF sample (Marshall 
1983), we estimate that the contribution of all QSOs brighter than B = 19.5 to the 
XRB is 21 ± 8%. This estimate, as well as that of point (a), relies only on data from 
complete samples, and therefore is not subject to the majority of uncertain assumptions, 
parameters, and extrapolations which are usually required in such an application. The 
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major uncertainty in these estimates is the value assumed for the X-ray spectral slope: 
we conservatively (if the estimates above are regarded as upper limits) adopt Qj = 0.5. 

(c) Because of their apparent dependence on optical luminosity and/or redshift, the pa- 
rameters < a ox > and aS?, as well as the distribution of a ox values, are directly inter- 
pretable, model independent quantities only if measured in a narrow region of L opl . z 
space. Our measurements for sample N (which is. by construction, such a “narrow’ 1 
sample) yield < a ox > = 1.58. and a% x = 1.50, for QSOs with < log (/ 2 5 oo) >= 31.44 
(Mb — —27.5) and < 2 >= 2.1. These values agree well with the predictions of 
Avni and Tananbaum (1982. 1986) for dependence of < a OI > on (primarily) optical 
luminosity. This conclusion, plus our model-independent, non-parametric, empirical 
comparison of the X-ray data of sample N with that of the recent ly published “XBQS’h 
do not provide support for the recent suggestion by Schmidt and Green (1986) that 
the primary dependence of a ox is on 2 . 

(d) The width of the a oz distribution for sample X, o(a Dr ) = 0.12 - 0.14, is smaller than 
that found by Avni and Tananbaum (1982) for a more heterogenous sample. Although 
this smaller width potentially could resolve some conflicts between inferences derived 
from X-ray versus optically selected QSOs, a discrepancy which in part has motivated 
Schmidt and Green (1986) to suggest reconsideration of the form of X-ray evolution of 
QSOs, it is likely that a ox threshold effects in survival analysis are artificially narrowing 
the width of the distribution. The half-width of the distribution below < a ox > is of 
order 0.1 for sample N, and thus, unless the actual distribution is quite skew, its true 
full width is probably also <7(a 0I ) ~ 0.2, in agreement with the Avni and Tananbaum 
value. 

(e) If we adopt recent QSO surface density counts and the hypothesis of an < a OI > 
dependence on L opi . we find that QSOs are capable of making the dominant contribu- 
tion to the XRB at soft energies. Contributions as large as 70 a are consistent with 
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these data, with the typical contributor a moderate redshift (2 ~ 1 — 2) QSO with 
B ~ 19 — 21. Considerable uncertainty is introduced into this analysis because of the 
poorly known distribution of X-ray spectral indices of QSOs. However, ~ 50% of the 
QSO contribution (regardless of whether or not it is the dominant one) comes from 
QSOs brighter than B ~ 20. 

We are indebted to the staffs at the Harvard-Smithsonian Center for Astrophysics, 
especially Drs. H. Tananbaum and F. Seward, and at KPXO, CTIO, and the CFHT for 
their aid in obtaining and reducing the data described here, and to a long list of Einstein 
Observatory guest observers who permitted us to examine portions of their data prior to 
public release. Dr. D. Helfand initially suggested the application of the stacking analysis 
to our problem. This work has been supported in part by NASA Grant NSG8-433; the 
first author also gratefully acknowledges the financial support of the ARCS Foundation. 


APPENDIX 

In many astronomical applications, the presence of detection thresholds results in a 
quite non-random distribution of the values of the non-detected quantities. We now show 
by Monte Carlo simulations that such ‘'threshold" effects and an assumption commonly 
made in survival analysis can produce some very misleading results. Note that thresholds 
in q oi units can correspond to joint thresholds in optical and X-ray flux sensitivities. 

Using a Gaussian random number generator, we generated a synthetic sample of 100 
objects which resulted in a distribution with < a OJ >= 1.55, 0 — 0.2, and ctoi = 1.45 
(A85); these values are similar to those found for sample N. The Gaussian form is not 
intended to necessarily represent the actual distribution, but is merely used as a convenient 
computational form. 

First, consider this distribution to be an idealized sample N distribution; in particular, 
suppose it to be the a ox distribution of 100 QSOs. all of the same L opt and 2 , and with all 
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objects imaged in X-rays to the same X-ray flux threshold. For this idealized sample, the 
X-ray flux limit directly translates into a uniform lower limit on a oz for the non-detections 
in the sample. We impose this uniform a oz threshold on the idealized distribution, and 
subject the resulting data to conventional survival analysis to determine < a 0I > and q| 5- 
Finally, we repeat this process for various limiting thresholds. Throughout, we follow the 
usual practice of treating the highest a ox observation (limit or detection) as though it were 
uncensored, an assumption shown below to be poor in certain circumstances. 

The results of this Monte Carlo simulation are shown in Table Al. The mean a ox value 
estimated from survival analysis is in quite serious disagreement with the true value for 
this extreme censoring (threshold) pattern, even when the threshold is sensitive enough to 
detect ~ 50% of the objects. The estimated from survival analysis is somewhat more 
reliable, due to its heavy weighting of X-ray bright objects. The systematic underestimate 
evident in Table Al can be understood as follows. Again following common practice, we 
have assumed that the survivor function S(a ox ) goes to zero at the highest observation 
(censored or uncensored). That is, the non-parainetric survival analysis artificially trun- 
cates the distribution; without some prior knowledge of the true distribution function, this 
is the best it can do. 

Note further from Table Al that such threshold censoring can change the apparent 
relation between a ox and a% x , and that the survival analysis error estimates are extremely 
small (and misleading!) in the presence of significant threshold censoring. This latter 
point is also readily understood in terms of the premature truncation of the estimated 
distribution: for a ox thresholds less than < a 0I >. there are very few objects in the parent 
distribution with a ox much below the threshold, and of course none (due to truncation) 
above that threshold, so that in the survival analysis estimated distribution, most objects 
fall close to the threshold. Thus the apparent width of the distribution is underestimated, 
and the distribution is artificially skewed toward low a OT values. 
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If a reliable guess cannot be made for the a ox value at which the true distribution has 
S(a 0I ) = 0, then the uncertainty in < a ox > in the positive direction becomes arbitrarily 
large in the presence of significant threshold censoring (Avni et al. 1980). Then, (see Table 
Al), in the presence of significant threshold censoring, it is more appropriate to consider 
survival analysis parameter estimates as lower limits on < a ox > and a-|f. 

Now we suppose instead that there are two idealized N-like QSO samples. Further, 
suppose that one sample, sample A, consists of objects at optical luminosity Z^qo and with 
redshift a. and the second sample, sample B, consists of objects with optical luminosity 
^2500 ar *d also at redshift 2 (assume that l^oo > ^ 2500 ) i ar *d suppose that both samples A 
and B have the same Gaussian parent distribution of a oz values as that described in the 
preceeding few paragraphs. By construction, we have assumed no dependence of a ox on 
L opt . Now we apply a uniform X-ray flux threshold to both samples; since sample A is a 
collection of more optically luminous objects, it has a more sensitive a ox threshold than 
sample B, i.e., ^ > ( ai™) . Conceptually, the results of Table Al also apply to the 

latter hypothetical situation; it is seen that the standard survival analysis assumption leads 
to the clearly mistaken conclusion that the < L x / L opi > ratio is higher in the lower optical 
luminosity sample B, than in the higher optical luminosity sample A. For a given redshift, 
only the optically (and hence X-ray) more luminous objects are detectable at sufficiently 
high a ox values that the shape of the actual distribution function is well sampled. 

Qualitatively, such spurious threshold effects could mimic the sort of dependence of a 0J 
on L op t given physical interpretation by, for example. Tucker (1983) and Netzer (1985). 
In practice, this sort of threshold effect competes against observational selection: most 
high L opt QSOs are also at high redshift, and apparently faint in the optical. This latter 
selection effect tends to make the a ox threshold for high L op t objects systematically smaller 
than for lower redshift. objects. Thus, one nmst consider the censoring pattern to assess the 
potential influence of threshold effects. We have done this (see A85) for the data analyzed 
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by Avni and Tananbaum (1982). and see no evidence that the dependence of a OI on L opt 
which they find results spuriously from the aforementioned threshold effect. Note also that 
the above sort of sharp cutoff threshold pattern arises naturally in other types of samples 
in addition to those like sample N (see A85). 

Finally, we note that an assumption fundamental to survival analysis is that detections 
and non-detections arise from the same parent population. This seems reasonable in X-ray 
studies of QSOs, since substantial fractions (e.g., 50% in the BF and Ivriss and Canizares 
samples, and > 80% in the XBQS) of optically selected QSOs are in fact X-ray sources. 
However, that this assumption has been made must be carefully kept in mind when using 
survival analysis. For example, were there two distinct populations of QSOs - one X- 
ray loud and the other X-ray quiet - survival analysis could lead to extremely erroneous 
conclusions about the distribution of X-ray properties of all QSOs, if it were assumed that 
X-ray quiet and X-ray loud QSOs alike came from the same parent distribution of X-ray 
sources: see A85 for a specific numerical example. Avni and Tananbaum (1986) have 
argued (in a somewhat distribution dependent fashion) that any populat ion of X-ray quiet 
QSOs is a small (< 10%) fraction of the entire QSO population. 

In summary, the effects discussed here, while potentially important in survival analysis, 
do not appear to have materially affected published results on the XRB. 

We have also performed Monte Carlo simulations to assess the influence of X-ray 
threshold effects on stacking. The model here simulates 100 QSOs, each imaged in X-rays 
for t, seconds (all t t 's are assumed equal). In imitation of sample N, the QSOs are assumed 
to all have z — 2.1 and log (/ 2 5 oo) = 31.44; a oz values are again distributed according to 
the Gaussian distribution used in the earlier Monte Carlo simulation. The ensemble mean 
IPC counting rate for this synthetic sample would be < $ >= 1.66 ± 0.096 x 10~ 3 counts 
s -3 in the 0. 5-3.0 keV band. For each QSO/image, noise is introduced into the individual 
S/’s. B 7 ~ s, and T t ' s. Then the 100 noise-added images are stacked and the ensemble values 
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of T and B for the stacked image are determined as described above. Finally, this whole 
process is repeated for different assumed s (i.e., different count rate thresholds) yielding 
the results of Table A2; the detection fraction listed there is the number of objects that 
would have been accepted as X-ray detections under guidelines given in § 3. The results 
are quite encouraging: even in the presence of substantial non-random censoring, stacking 
apparently is not strongly affected by threshold effects. Compare Table A2 for stacking, 
with Table A1 for survival analysis, and note that an approximate conversion between 
an a ox threshold and an integration time threshold can be inferred by comparing the 
“detection fraction” columns in the two tables. 
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Table 1. The Complete QSO Sample ‘•C■ , ( B< 19.5, 1.8 < z < 3.0) 


Object 0 

a (1950) 

<5 (1950) 

B 

log (^ 2500 ) 

z 6 

08366+654 

08 /, 36 m 37.5* 

+65°24'0T' 

19.0 

31.57 

1.9 

08388+131 

08 38 48.4 

+ 13 10 00 

19.1 

31.53 

1.88 

08388+133 

08 38 50.5 

+ 13 19 42 

19.0 

31.52 

1.80 

08500+283 

08 50 03.6 

+ 28 18 32 

19.3 

32.11 

2.90 

09038 + 167 

09 03 52.6 

+ 16 46 15 

18.2 

32.14 

2.4 

09379+121 

09 37 56.8 

+ 12 09 33 

19.0 

31.98 

2.7 C 

09388+117 

09 38 50.8 

+ 11 42 54 

18.6 

31.87 

2.2 

09392 + 121 

09 39 14.5 

+ 12 07 24 

19.0 

31.46 

1.80 

11131 + 183 

11 13 08.4 

+ 18 21 18 

18.6 

31.83 

2.20 

11156+180 

11 15 33.3 

+ 18 02 39 

18.1 

31.85 

1.90 c 

12071 + 398 rf 

12 07 11.6 

+ 39 53 24 

19.4 

31.61 

2.33 c 

12076+399 

12 07 38.6 

+ 39 56 1 i 

17.5 

32.39 

2.4 

12302+120 

12 30 11.0 

+ 12 02 59 

19.1 

31.45 

1.90 

12336+264 

12 33 28.7 

+ 26 29 54 

19.1 

31.75 

2.40 

12354 + 264 rf 

12 35 25.8 

+26 27 30 

19.4 

31.68 

2.5 

14149+251 

14 14 59.5 

+25 06 02 

18.7 

31.63 

1.87 

14151 + 254 

14 15 08.3 

+ 25 28 45 

19.3 

31.64 

2.31 

15181 + 201 

15 18 07.6 

+20 07 50 

19.2 

31.62 

2.1 

16017+184 d 

16 01 47.0 

+18 25 34 

19.1 

31.51 

1.94 c 

17272+502 

17 27 12.5 

+50 15 07 

19.1 

31.67 

2.1 

17272+499 

17 27 17.0 

+49 56 41 

19.3 

31.44 

1.9 

17274 + 503 

17 27 28.1 

+50 21 03 

19.2 

31.66 

2.2 

20382-012 

20 38 16.6 

-01 16 21 

19.1 

32.21 

2.9 

21342-149 

21 34 13.1 

-14 55 56 

IS. 3 

32.00 

2.20 
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21357-147 

21 35 45.7 

-14 45 10 

19.3 

31.57 

2.1 

21556-302 

21 55 40.4 

-30 14 27 

18.9 

31.75 

2.20 

21570-302 

21 57 05.9 

-30 14 52 

18.4 

31.77 

1.9 

Obj. 5 



18.3 

31.91 

1.96' 

Obj. 6 



19.3 

31.67 

2.20 

Obj. 7 



19.2 

31.66 

2.2 

Obj. 8 



19.1 

32.23 

3.00 

Obj. 10 



19.3 

31.53 

2.06 

Obj. 11 



19.1 

31.47 

1.90 

Obj. 12 



19.2 

31.41 

1.90 

Obj. 15 



19.1 

31.60 

2.0 

Obj. 16 



18.2 

32.32 

2.7 

Obj. 17 



19.0 

31.93 

2.6 

Obj. 19 



19.4 

31.63 

2.10 

Obj. 21 



19.4 

31.57 

2.10 


a The nomenclature “Obj." is used here and in subsequent tables to honor the propri- 
etary data rights of other observers who generously loaned to us certain fields for this 
investigation. 

b Redshifts here and in Table 2 are based on identification of strongest feature as La: 
when more than a single line on the grism/grens plate is detected at a consistent 
redshift (usually C IV), two decimal places are listed for the redshift (although redshift 
is still uncertain to ~ 0.1 in z). 

c Slit spectrum confirms redshift . 


d Previously cataloged objects (see references in Veron-Cetty and Veron 19S5): 12071 -f 
398 = 1207 + 3911 "3, 12354 + 264 = Wee 73. 16017 + 184 = 1601 4- 184. 



Table 2. The “Narrow" QSO Sample “N” 
(31.0 < log (/ 2500 ) < 32.0, 1.8 < 2 < 2.5) 


Object 

a (1950) 

6 (1950) 

B 

l°g (^ 2500 ) 

z 

00151 + 160 

00 b 15 m 09.9 s 

+ 16°03'13" 

19.9 

31.33 

2.20 

00159 + 155 

00 15 54.8 

4-15 35 48 

20.6 

31.07 

2.30 

01121 - 014 

01 12 06.6 

-01 26 45 

20.3 

31.17 

2.20 

0130 - 404° 

01 30 45.9 

-40 25 49 

19.1 

31.69 

2.16 

0130 - 406° 

01 30 50.7 

-40 3S 13 

18.9 

31.88 

2.39 

0131 - 405° 

01 31 09.0 

-40 35 44 

20.4 

31.21 

2.25 

0131 - 402 a 

01 31 36.6 

-40 17 56 

20.2 

31.20 

2.11 

0131 - 409° 

01 31 40.0 

-40 54 23 

19.7 

31 .55 

2.36 

0132 - 403 a 

01 32 09.6 

-40 18 08 

19.4 

31.58 

2.18 

0132 - 408 a 

01 32 20.8 

-40 48 26 

20.7 

31.23 

2.46 

0132 - 406° 

01 32 22.2 

-40 37 18 

19.4 

31.56 

2.15 

02031 + 151 

02 03 07.3 

4" 15 09 0 1 

20.3 

31.10 

2.00 

02031 + 152 

02 03 08.4 

4-15 17 04 

19.7 

31.50 

2.38 

02036 + 150 

02 03 38.4 

+ 15 01 05 

19.9 

31.32 

2.10 

02057 + 150 

02 05 44.8 

+ 15 00 41 

20.5 

31.20 

2.40 

03057 + 172 

03 05 45.9 

+ 17 12 12 

20.4 

31.27 

2.40 

03061 + 169 

03 06 06.7 

+ 16 54 23 

20.3 

31.21 

2.14 

03074 + 172 

03 07 28.9 

+ 17 16 37 

19.2 

31.66 

2.28 

0S366 + 654 b 

08 36 37.5 

+65 24 01 

19.0 

31.57 

1.9 

08388 + 131 6 

08 38 48.4 

+ 13 10 00 

19.1 

31.53 

1.8S 

08388 + 133 b 

08 38 50.5 

+ 13 19 42 

19.0 

31.52 

1.80 

09382 + 117 

09 38 13.9 

+ 11 45 32 

19.9 

31.37 

2.30 



09382 + 120 


09 38 14.7 +12 

09383 + 120 09 38 23.9 +12 

09388 + 117 6 09 38 50.8 +11 

09392 + 121 b 09 39 14.5 +12 

09392 + 117 09 39 17.2 +11 

11131 + 183 6 11 13 08.4 +18 

11136 + 182 11 13 41.9 +18 

11143 + 184 11 14 19.2 +18 

11156 + ISO 6 11 15 33.3 +18 

11147 + 183 11 14 46.6 +18 

12071 + 398 b 12 07 11.6 +39 

12275 + 024 12 27 35.3 +02 

12292 + 116 ° 12 29 12.4 +11 

12302 + 120 b 12 30 11.0 +12 

12336 + 264 6 12 33 28.7 +26 

12354 + 264 b 12 35 25.8 +26 

14144 + 256 14 14 29.9 +25 

14148 + 252 14 14 49.7 +25 

14149 + 251 b 14 14 59.5 +25 

14151 + 254 b 14 15 08.3 +25 

15090 - 092 15 09 03.4 -09 

15100 - 089 15 10 00.5 -08 

15181 + 201 b 15 18 07.6 +20 

15184 + 202 15 18 25.9 +20 

16016 + 184 16 01 38.3 +18 

16017 + 184 b 16 01 47.0 +18 

16 23 25.4 +26 
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05 45 

19.5 

31.42 

2.00 

05 58 

20.2 

31.32 

2.4 

42 54 

18.6 

31.87 

2.2 

07 24 

19.0 

31.46 

1.80 

46 02 

19.8 

31.21 

1.9 

21 18 

18.6 

31.83 

2.20 

12 53 

19.5 

31.29 

1.9 

28 38 

20.3 

31.15 

2.20 

02 39 

18.1 

31.85 

1.90 

19 42 

19.7 

31.21 

1.9 

53 24 

19.4 

31.61 

2.33 

28 47 

20.2 

31.12 

2.0 

37 48 

19.7 

31.14 

1.80 

02 59 

19.1 

31.45 

1.90 

29 54 

19.1 

31 .75 

2.40 

27 30 

19.4 

31.68 

2.5 

36 42 

19.6 

31.22 

1.80 

13 34 

19.7 

31.20 

1.83 

06 02 

18.7 

31.63 

1 . S7 

28 45 

19.3 

31.64 

2.31 

17 09 

19.9 

31.18 

1.9 

57 48 

20.4 

31.17 

2.1 

07 50 

19.2 

31.62 

2.1 

17 33 

19.9 

31.34 

2.10 

24 08 

20.2 

31.24 

2.31 

25 34 

19.1 

31.51 

1.94 

54 41 

19.6 

31.66 

2.44 


1623 + 269 ° 
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1624 + 269° 

16 24 03.1 

+ 26 57 33 

18.3 

32.00 

2.20 

17252 + 499 

17 25 14.1 

+49 57 46 

19.9 

31.20 

1.90 

17257 + 503 

17 25 43.4 

+50 18 15 

20.4 

31.15 

2.1 

17264 + 504 

17 26 25.3 

+50 28 30 

19.9 

31.20 

1.9 

17272 + 502 fc 

17 27 12.5 

+ 50 15 07 

19.1 

31.67 

2.1 

17272 + 499 6 

17 27 17.0 

+49 56 41 

19.3 

31.44 

1.9 

17274 + 503 fc 

17 27 28.1 

+ 50 21 03 

19.2 

31.66 

2.2 

17449 + 206 

17 44 59.1 

+ 20 36 05 

19.0 

31.91 

2.41 

17465 + 201 

17 46 33.3 

+20 10 58 

19.0 

31.64 

1.90 

20374 - 007 

20 37 25.8 

-00 44 39 

19.9 

31.47 

2.3 

20381 - Oil 

20 38 10.9 

-01 09 10 

20.4 

31.25 

2.26 

21189 + 168 

21 18 59.9 

+ 16 51 26 

19.5 

31.64 

2.30 

21259 - 148 

21 25 57.4 

-14 51 02 

19.9 

31.38 

2.30 

21265 - 150 

21 26 34.1 

-15 02 39 

20.3 

31.18 

2.1 

21342 - 149 fc 

21 34 13.1 

— 14 55 56 

18.3 

32.00 

2.20 

21354 - 145 

21 35 27.8 

-14 32 11 

19.9 

31.18 

1.90 

21357 - 147 6 

21 35 45.7 

-14 45 10 

19.3 

31.57 

2.1 

21416 + 037 d 

21 41 40.5 

+03 45 0/ 

20.2 

31.05 

1.8 

21431 + 040 

21 43 10.8 

+ 04 03 59 

19.2 

31.57 

2.0 

21558 + 034 

21 55 51.2 

+03 29 33 

19.0 

31.56 

1.9 

21556 - 302 6 

21 55 40.4 

-30 14 27 

18.9 

31.75 

2.20 

21570 - 302 6 

21 57 05.9 

-30 14 52 

18.4 

31.77 

1.9 

Obj. 1 



20.4 

31.13 

2.20 

Obj. 2 



18.7 

31.97 

2.5 

Obj. 3 



20.3 

31.20 

2.30 

Obj. 4 



19.1 

31.61 

2.0 


Obj. 5 b 
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Obj. 

6 6 




Obj. 

7 b 




Obj. 

9 

20.2 

31.18 

2.1 

Obj. 

10 6 




Obj. 

ll 6 




Obj. 

12 b 




Obj. 

13 

20.1 

31.22 

2.1 

Obj. 

14 

20.2 

31.20 

2.2 

Obj. 

15 b 




Obj. 

18 

18.9 

31.81 

2.0 

Obj. 

19 b 




Obj. 

20 

20.2 

31.35 

2.3 

Obj. 

2l 6 





a QSOs discovered in previous grism surveys, whose X-ray properties have been studied 
by Kriss and Canizares (1985): l 2 500 has been recalculated assuming o c = 1.0. 

b QSO also in Sample C: details in Table 1. 

c Identical to object U9 of He et al. (1984). 

d Coincident within the positional errors with the 5 GHz radio source MG2144+0358 
(Bennett et al. 1986). 
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Table 3. X-Ray Data for Sample C and N QSOs° 


Fobs X 10 13 

(erg cm -2 s _1 , 


Object 

Significance 6 

0.5-3.0keV) 

log O 2 ke\ ; ) 

&OI 

00151 + 160 

0.7 

< 0.688 

< 28.24 

> 1.19 

00159 + 155 

0.2 

< 0.565 

< 27.70 

> 1.29 

<j 

o 

1 

o 

CO 

1—1 

o 



27.97 

1.43 

0130 - 406 c 



< 27.62 

> 1.64 

0131 - 405 c 



< 27.54 

> 1.41 

0131 - 402 c 



27.77 

1.32 

0131 - 409 c 



< 27.58 

> 1.52 

0132 - 403 c 



27.83 

1.44 

0132 - 408 f 



< 27.56 

> 1.41 

0132 - 406 c 



< 27.44 

> 1.58 

02031 + 151 

-0.7 

< 0.641 

< 27.60 

> 1.34 

02036 + 150 

1.3 

< 0.633 

< 27.65 

> 1.41 

03061 + 169 

2.1 

< 0.645 

< 27.71 

> 1.34 

08366 + 654 

0.4 

< 0.728 

< 27.59 

> 1.53 

08388 + 131 

2.1 

< 0.766 

< 27.60 

> 1.51 

08388 + 133 

1.7 

< 0.735 

< 27.53 

> 1.53 

08500 + 283 

0.2 

< 0.692 

< 28.07 

> 1.55 

09038 + 167 

0.0 

< 0.882 

< 27.95 

> 1.61 

09379 + 121 

-0.2 

< 0.601 

< 27.91 

> 1.56 

09382 + 117 

-0.1 

< 0.670 

< 27.77 

> 1.38 

09382 + 120 

1.3 

< 0.576 

< 27.54 

> 1.49 

09383 + 120 

-0.8 

< 0.597 

< 27.77 

> 1.36 
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+ 147 

0 9392 + 121 

09392 + 

477*7 + 7* 

\w& + 182 

47743 + 7 * 

U 147 + 183 

\\\& + l8 ° 
1207 ^ 398 
yZO'O A 3 " 
tfl73 + 034 
12 292 + 

}O302 + 72° 
yo3 ' o4 + 264 
X 414S + 232 
!4149 + 251 


2-0 
4.0 
2.2 
-o.e 
- 2-0 
2-7 
4. a 
4-9 
4-S 
0A 
-0.0 
2-0 
-OA 
4A 
2-4 
-0.9 
4-0 


^ 0.080 
< 0-020 
< 0-002 
< 0-090 
< 0-082 
< 0-007 
< 0-092 
< 0-790 
< 0.008 


< 27-73 
< 27-40 
< 37 0 4 
< 27-00 
< 27 04 
<27-7° 

< 37 .00 
<37.07 


< 0-080 
^ 3.690 
< 0-747 
< 4.370 

< 0 -099 
< 0 -002 
< 4-700 
< 0-003 
< 2 - 07 ° 


> 4 -09 
> 4-04 
> 7-47 
> 4-00 
> 4-44 
> 4.33 
>4-40 
> 4-03 
>7-49 
>7.70 
> 7 A 2 
>4.39 
>4.39 


< 27.73 
< 27 -87 
<38.20 
<37 03 
< 27.83 
< 27 -87 
< 27 -47 
< 29 -49 
< 27 -80 
<28-70 


14151+ 254 

4-0 

< 4-200 

1&090 ' ° 9 ' 

0-2 

< 0-072 

75700 - 089 

- 0-9 

< 0-620 

15 ^ + 2 ° 2 

3.3 

< 0-073 

16016 + 184 

4-0 


I 60 l " i +l84 



1633 + 269 


< 0-094 

1624 + 269 

- 2.0 

< 0-008 

+ 409 

\ 72 o 2 ^ 

- 0-3 

< 0-004 

, - l 303 
\720 * + 

4-4 



< 28.37 
< 27 00 
<27 -70 
< 27 07 
< 27 03 

37 .0° 

< 27 -00 
< 27 -09 
< 27 A8 


> 4-49 
> 4-40 
> 0.82 
> 4-40 
> 4-70 
> 4-07 
> 4-40 
>4-34 
> 4-04 
> 4-00 
4-70 
> 4-42 
>7.37 
> 4-43 


17264 + 504 



17272 + °°2 


< 87.7 

— 1.2 


4S 

< 1.53 X 10 4 

< 0.633 

< 0.521 


< 30.03 

< 27.53 
<27.61 


> 0-63 

> 1.50 

> 1.55 

> 1.18 


X7274 + 503 

l\ 

0 ^ 

w 

< 0-646 

< 20 .3° 

> 1-34 

< 27. <6 

17465 + 201 

1.5 

< 0.643 

> 1.62 

< 28.00 > 

20381 " 011 

2.1 

< 0.561 

> 1-38 

< 27.79 >A 

20382 - 012 

1.6 

< 0.674 

- R7 > 135 

< 2/ -6' 

21259 - 148 

0.6 

< 0.662 

> 1 63 

<27.o 

21265 - 1°° 

1.3 

< 0.710 

> 1.36 

<27.63 

21342 - l 49 

1.2 

< 0.787 

> 1.49 

< 27.10 

21354 - I 45 

2.0 

< 0.715 

> 1 36 

< 2< -o2 

21357 - 1 4< 

1.2 

< 0.658 

<27.72 > 145 

21416 + 03 ‘ 

-1.6 

< 0-836 

< 27.59 > l '°' 

21431 + 040 

0.7 

< 0.712 

> 1.46 

< 27-94 

21558 + 034 

< 2-1 

< 1.130 

> 1-54 

<27.0 

21556 - 302 

1.7 

< 1.080 

- -a 136 

2< .o9 

21570 - 302 

2-7 

0.493 

1.56 

27.90 

Obj. 1 

3.0 

0.689 

1.20 

28.06 

• > 
& 
0 

6.2 

1.280 

1.54 

27.59 

Obj. 3 

3.6 

0.597 

1.65 

27.61 

Obj. 4 

4.8 

0.533 

1.60 

27.49 


Obj. 6 
Obj. 6 
Obj. 7 
Obj- 8 
Obj. 0 
Obj. 1° 
Obj. H 


2.8 

3.4 

4.4 
2.9 
3.5 
2.8 


0.925 

1.110 

0.685 

0.765 

0.265 


27.86 

28.33 

27.68 

27.69 
27.60 


1.46 

1.50 

1.35 

1.48 

1.48 
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Obj. 

12 

3.3 

0.876 

27.66 

1.44 

Obj. 

13 

3.4 

0.914 

27.79 

1.32 

Obj. 

14 

4.5 

1.220 

27.97 

1.24 

Obj. 

15 

2.1 

< 0.648 

< 28.06 

> 1.36 

Obj. 

16 

-0.2 

< 0.459 

< 28.27 

> 1.56 

Obj. 

17 

2.7 

0.664 

27.92 

1.54 

Obj. 

18 

10.1 

2.050 

’ 28.13 

1.41 

Obj. 

19 

2.9 

0.693 

27.70 

1.51 

Obj. 

20 

7.9 

2.800 

28.41 

1.13 

Obj. 

21 

2.7 

0.688 

27.74 

1.48 


a QSOs in the sample which lie within an optical grism/grens plate, but lie outside the 
usable area of the X-ray image ( e.g on the “ribs" in the IPC). have incomplete X-ray 
information and are omitted from the table. 

b Ratio of source counts/(background counts) 1 / 2 in 3' radius circle. 

c X-ray data from Kriss and Canizares (1985). but with a OI adjusted for the assumption 
of a 0 = 1.0. 
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Table 4. Fractional Contribution of QSOs to the XRB, 
Assuming a ox Dependence on L opt 


Zeff — 
(B ma x) = 

0.5 

(19.87) 

1.0 

(21-77) 

1.5 

(22.99) 

2.0 

(23.90) 

2.5 

(24.64) 

3.0 

(25.27) 

3.5 

(25.81) 

A B 








16-17 

0.006 

0.004 


0.003 


0.002 


17-18 

0.025 

0.017 


0.012 


0.009 


18-19 

0.109 

0.074 


0.049 


0.039 


19-20 

0.363 

0.316 


0.212 


0.166 


20-21 


0.206 


0.138 


0.108 


21-22 


0.106 


0.091 


0.071 


22-23 




0.078 


0.061 


23-24 




0.060 


0.052 


24-25 






0.045 


25-26 






0.011 


16-B max 

0.503 

0.722 

0.686 

0.642 

0.600 

0.564 

0.532 

16-19.5 

0.293 

0.198 

0.157 

0.132 

0.116 

0.104 

0.095 



Table Al. Non-Random Threshold Effects in Survival Analysis 


detection 


,-thresh 

a ox 

"C cx OT > 

n .eff 

fraction 

1.3 

1.29 ± 0.002 (13.0) 

1.29 ±0.01 (7.2) 

.11 

1.4 

1.38 ± 0.005 ( 8.2) 

1.37 (—.01, ±.005) (3.6) 

.27 

1.5 

1.44 ± 0.009 ( 5.0) 

1.41(-0.01. ±.02) (1.4) 

.44 

1.6 

1.49 ± 0.013 ( 2.5) 

1.44( — .02. ±.01) (< 1) 

.62 

1.7 

1.52 ± 0.016 ( 1.2) 

1.44( — .01. ±.02) (< 1) 

.78 

1.8 

1.54 ± 0.018 (< 1) 

1.45 ± 0.02 (< 1) 

.90 

1.9 

1.54 ±0.019 (< 1) 

1.45 ± 0.02 (< 1) 

.93 

2.0 

1.55 ± 0.020 (< 1) 

1.45 ± 0.02 (< 1) 

.99 

2.1 

1.55 ± 0.020 (< 1) 

1.45 ± 0.02 (< 1) 

1.00 


Note: numbers in the rightmost set of parentheses within the second and third columns 
show the significance (in o' s) of the difference between the survival analysis estimate and 
the chosen values of < a oz > = 1.55 and = 1.45. 
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Table A2. Threshold Effects in Stacking 


t«(s) 

< s > (10 3 counts s ! ) 

detection 

fraction 

20,000 

1.634 ±0.090 

0.38 

10,000 

1.624 ±0.128 

0.27 

5,000 

1.610 ±0.179 

0.17 


1.000 


1.551 ± 0.396 


0.03 
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Figure Captions 

Figure 1. The distribution of apparent B magnitudes for the grism/grens emission line 
selected QSOs (mainly at high redshift) discovered on the best quality plate material. Note 
the decrease at magnitudes fainter than the peak, which we interpret as due to the onset 
of incompleteness at B ~ 19.5. 

Figure 2. The reliability of magnitudes estimated from POSS has been tested through 
measurement of two samples with photometric data. Triangles , the spectrophotometrically 
calibrated QSO sample of Vaucher and Weedman (1980). using those authors’ 4500 A 
continuum magnitudes; squares, the M15 photoelectric sequence of Sandage and Katem 
(1977). 

Figure 3. The survivor function for a OI values of QSOs in sample N (1.8 < z < 2.5. 31.0 < 
/ o<7 (/2500 ) — 32.0). This function S(a) describes the probability that a QSO will have an 
a 0I value in excess of a. Vertical broken lines, values typical of X-ray selected samples 
(Gioia et al. 1984. Margon, Downes, and Chanan 1985). and the UY-selected sample 
studied by Marshall et al. (1984). Note that many of the QSOs studied here have higher 
a 0T values, and are thus significantly more "X-ray quiet", than those of past work. For 
our sample. aS? = 1.50 ± 0.03. 

Figure 4. A comparison of the sample N a PT survivor function (solid line) for high redshift 
QSOs with those of low redshift. XBQS QSOs of comparable L opt . Sample N QSOs 
( z — 1.8 — 2.5) span a range in optical luminosity of 30.8 < log(l 2500 ) < 32.0 and 30.5 < 
log (1 9500 ) < 31.6 for qo — 0 and 0.5. respectively. This function differs from that shown in 
Figure 3 because here we assume o c = 0.5 for ease of comparison with the XBQS. Broken 
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line. a ox distribution for the XBQS subset, that spans the same L opt range as sample N 
assuming q o = 0, and that occupies a similarly narrow, but much lower, redshift range. 
z = 0.16 — 0.62: this first XBQS subset includes 12 objects, of which 10 are detected in 
X-rays. Dotted line, a ox distribution of the XBQS subset in the appropriate L opt range 
for comparison with sample N, for the case of qo = 0.5: this subset spans 2 = 0.14 — 0.62, 
and includes 19 objects, 17 of which are detected in X-rays. A typical error bar for the 
XBQS distribution is shown in the upper right hand corner. The similarity of the survivor 
functions for sample N at < z >= 2.1 and the XBQS subsets at < z >— 0.3 — 0.4 (see text 
for statistics) argues against substantial redshift evolution of a ox . 

Figure 5. Survivor function for a ox values of sample A, with corrections applied for the 
correlation of a ox with optical luminosity. As described in the text, the corrections are 
sufficiently small that this distribution is virtually indistinguishable from that of Figure 3. 
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